Mesoscopic and microscopic dipole clusters: Structure and phase transitions 
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Two dimensional (2D) classical system of dipole particles confined by a quadratic potential is stud- 
ied. This system can be used as a model for rare electrons in semiconductor structures near a 
metal electrode, indirect excitons in coupled quantum dots etc. For clusters of A'' < 80 particles 
ground state configurations and appropriate eigenfrequencies and eigenvectors for the normal modes 
are found. Monte Carlo and molecular dynamic methods are used to study in detail the order - 
disorder transition (the "melting" of clusters). In mesoscopic clusters {N < 37) there is a hierarchy 
of transitions: at lower temperatures an intershell orientational disordering of pairs of shells takes 
place; at higher temperatures the intershell diffusion sets in and the shell structure disappears. In 
"macroscopic" clusters (A'' > 37) an orientational "melting" of only the outer shell is possible. The 
most stable clusters (having both maximal lowest nonzero eigenfrequencies and maximal temper- 
atures of total melting) are that of completed crystal shells which are concentric groups of nodes 
of 2D hexagonal lattice with a number of nodes placed in the center of them. The picture of dis- 
orderung in clusters is compared with that in infinite 2D dipole system. The study of the radial 
diffusion constant, the structure factor, the local minima distribution and other quantities shows 
that the melting temperature is a nonmonotonic function of the number of particles in the system. 
The dynamical equilibrium between "solidlike" and "orientationally disordered" forms of clusters is 
considered. 

PACS numbers: 61.46. +w, 68.65.+g, 36.40.Ei 
I. INTRODUCTION 



Properties of a set of physical systems can, under cer- 
tain conditions, be studied with the help of the model 
system of interacting dipoles. Particularly, it is a dipole 
- dipole interaction that is of main importance wJien 
small dielectric particles on a|-surface of electrolyteEl or 
monolayers of adsorbed atomscl are considered. Another 
class of the systems in which the most part of inter- 
esting properties is due to a dipole character of inter- 
action is the class of electron systems in semiconduc- 
tors of small size located in the vicinity of a metal elec- 
trode electrostatic image forces can play an essential role. 
This role can become crutial, for example, for semicon- 
ductors with a small width of the forbidden band, for 
which ima|Se forces can induce a semiconductor - metal 
transition.Ll Indirect excitons with spatially separated 
electrones and holes in coupled quantum wells also in- 
teract by a dipole law, leaduie to the crystallization at 
intermediate concentrations .qEI 

Properties of infinite two dimensional dipole sjtstems 
have been intensively studied both tieoreticallyQ and 
with the help of computer simulationglll (see also Ref. [ 
P). But the study of mesoscopic systems of interacting 
dipoles (dipole clusters) is also of a particular interest. 
Recent advances in the microlithography and one - elec- 
tronics revived markedly an interest to systems of small 
number of particles, interesting with the strong struc- 
tural sensitivity to the number of particles and unusual 
rearrangements with increasing temperatur£. Particu- 
larly in recent experiments of Ashoori et aLcl the charg- 
ing of a quantum dot was studied by the single electron 
capacitance spectroscopy method. One of the most in- 
triguing results of these experiments was the observation 
of an unusual dependence of the chemical potential on 



the particle number N , when the total energy of the dot 
as a function of N had a quasi - periodic component of 
a universal shape (i.e. the electron addition are group 
into bunches). This dependence was shown to be caused 
by geometric effects associated with a packing-of the tri- 
angular lattice into the circular confinement .113 It is ob- 
vious that if the dot of a large enough size is situated 
near a single metal gate the electrostatic image forces 
will play an important role and the long-range electron 
- electron Coulomb interaction potential will transform 
to the short-ranged dipole - dipole one between "com- 
posite" particles "electron -I- image charge". A plane 
microtrap of cooled atoms with (normal to the plane) 
dipole momenta induced by static electric field or by an 
electromagnetic standing wave could be another physical 
example of a dipole cluster. The value of the dipole mo- 
menta of trapped atoms can be controlled by varying the 
magnitude of the field or by tuning the field frequency. 
Yet another interesting example of a dipole cluster is a 
system of indirect excitons in vertically coupled semicon- 



ductor quantum dots (see e.g. Refs. [ 11, ij]). Under 
certain conditions the system of indirect excitons can be 
considered as a classical dipole cluster in the parabolic 
confinement (see below). Such double dots systems are 
within the scope of current fabrication technology. 

The consideration of different properties of dipole. 
Coulomb and logarithmic clusters provides also a possi- 
bility to study the role of the interaction potential in the 
shell structure of ground state configurations (at T = 0) 
and in the hierarchy of melting temperatures. One of 
the most interesting tasks is the search for "magic" , most 
stable clusters. By now the structure of magic number 
clusters has been most extensively studied by the exam- 
ple of Lennard - Jones [LJ) systepis pUjiat correspond. 



e.g., to noble gases atomic clusters.El EJ Of some inter 



1 



est is the question: does the set of distinctive properties 
of " magic" clusters correlate with the particular form of 
an interaction potential? It is important to realize the 
manner in which changing the type of interaction poten- 
tial from two - parameter in L J system to repelling one - 
parameter in a dipole system may reflect on the proper- 
ties of appropriate "magic" clusters, on the character of 
disordering effects and on the possibility of coexistence 
phenomena. 

In this paper we study in detail the ground state struc- 
ture of 2D dipole clusters consisting of a finite number 
of particles and the picture of their melting with increas- 
ing temperature. It turnes out that on the basis of the 
character of disordering ("melting") phenomena dipole 
clusters can be divided into "mesoscopic" and "macro- 
scopic" groups. In "mesoscopic" clusters of small number 
of particles {N < 37) there are two stages of disorder- 
ing: orientational intershell disordering of different pairs 
of shells and, at greater temperatures, radial disordering 
when particles begin to interchange between shells. The 
spectral analysis of ground state configurations of small 
clusters shows that the value of the minimal frequency 
of normal vibrations correlates with characteristic tem- 
peratures of orientational disordering. Increasing of the 
temperature of large particle number clusters {N > 37) 
leads to radial disordering and to disruption of the shell 
structure at a temperature which is a function of the 
distance from the center of the system. An analysis of 
a behavior of different thermodynamic properties, when 
the nonuniformity of the density is taken into account, 
enables us to make a conclusion about the strong de- 
pendence of the total melting temperature on a cluster 
position in a periodic Mendeleev - type table (i.e. on 
the distribution of particles over crystal shells of differ- 
ent symmetry). As an intermediate between mesoscopic 
and macroscopic systems the cluster of 37 particles is 
considered. The most interesting feature of this cluster 
is the presence of the temperature region of coexistence 
of "solidlike" and orientationally disordered forms. This 
behavior of the system with increasing temperature is 
much alike the picture of structural transitions in magic 
number atomic clusters. 

The paper is organized as follows. In Sec. |l| we de- 
scribe the model and briefly outline the methods that 
where used to find global minima configurations and to 
calculate different thermodynamic properties. Sec. [II 



is devoted to the description and discussion of ground 
state configurations. We consider the "faceting" of clus- 
ters, the eigenvectors and eigenfrequencies for the normal 
modes. In Sec. IV we present the results from Monte 
Carlo (MC) and molecular dynamic (MD) simulations at 
different temperatures and system sizes. We consider two 
cases: "mesoscopic" {N < 37, se e Sec. |IV A ) and "macro- 
scopic" {N > 37, see Sec. IV B) clu sters. Separately, as 
an intermediate case, in Sec. IV C we study the cluster 
of iV = 37 particles. Our conclusions are presented in 
Sec. 0. 



II. THE MODEL. NUMERICAL SIMULATION 

The exposition below concerns a cluster of indirect ex- 
citons in double semiconductor dots, but the analysis can 
be obviously extended to any other above - mentioned 
dipole system. Let us consider two vertically coupled 2D 
semiconductor dots of a characteristic linear size L that 
are spaced at h apart. Under laser pulse on the forbidden 
band frequency, a formation of indirect excitons that are 
bounded states of electrons and holes belonging to dif- 
ferent dots is possible (see. works [ 0,|2| and references 
therein). In this work we are interested in the case that 
the exciton radius a{h) is much smaller of the mean dis- 



tance LjyN between them; a(/i) ~ Oq at the interdot 
separation h < Gq and a{h) ~ Oq^^^/i'^/^ at ft, ^ where 
Oq — K / (Am* e^) is effective Bohr radius of 2D exciton 
in a media with dielectric constant k, N stands for the 
(mean) number of excitons in the system. In addition 
we suppose that quantum effects in the cluster can be 
neglected in view of the condition k^T ^ N / {m* L"^) . 

Provided that conditions discussed above are fulfilled, 
excitons in quantum dots can be considered as 2D classi- 
cal system of particles with unidirected dipole momenta 
d ~ 2eh/K. A confinement potential can be chosen in 
a parabolic form with strength a. The Hamiltonian for 
such a system have the form 



H = 



r, - r. 



N 

1=1 



The Hamiltonian can be written in a dimensionless form 
if we express the coordinates and energy in the following 



units: ro = d^l^ ja^l^, Eo 
Hamiltonian becomes 



H 



1 



N 

P+Ei 



In such units the 



(1) 



From here on all the results will be given in the units in- 
troduced above. When considering thermodynamic prop- 
erties of the system we will use the dimensionless tem- 
perature T = khT/ar^. 

To find out ground state configurations of (|^) we ap- 
plied the following methods: j— , 

1) Classical simulated annealing method (CSA).Ej 

2) Combined Monte Carlo + gradient search" method. 

All of the results presented below (see Table I) were in- 
dependently obtained with the help of these algorithms in 
order to improve the reliability of results. Of course, none 
of existing algorithms of searching for the global mini- 
mum of multidimensional function can guarantee that a 
configuration obtained is the global minimum one. In or- 
der to overcome this difficulty we considered as many as 
200 random initial configurations. This approach made 
it possible to study local minima as well as appropriate 
regions of catchment ("relative weights" of local minima). 
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In view of intrinsic statistical nature of the CSA 
method, the problem of localization of the system in lo- 
cal minima is solved in this approach much more easy 
than in different gradient methods. When using CSA, 
the system is modelled at some artificially introduced 
temperature T{t) which is gradually decreased with the 
time t of experiment. This make it possible to sim- 
ulate the thermal noise and delocalize a system from 
metastable states. By starting from sufficiently high 
temperature of T{t = 0) = 10, at each t we per- 
formed ^200 MC steps per particle, whereupon the tem- 
perature was decreased as T{t + 1) = 0.98T(t). The 
trial move x ^ x + dx ai each MC step was gener- 
ated in accordance with Gaussian probability distribu- 
tion Il{Sx) ^ exp {—Sx'^ /T(t)). This move was accepted 
with the probability p = min{l, exp (- A£'/T(i))}. To 
gain in the accuracy of results several hundreds of con- 
jugate gradient steps were performed at the last stage of 
the algorithm. 

It turned out that CSA is not as fast the algorithm as 
the combined MC + gradient search method. A search 
for minima within the last approach consists of repeated 
(up to 10** combined moves) applying of the following 
methods: ■— ■ 

1) Classical random searchE3 when the trial move 
X ^ X + 5x IS accepted if it reduces the energy. The 
maximal size of the trial move max{|5a;|} was chosen au- 
tomatically to assure the acceptance probability of a new 
configuration of 0.1. 

2) Gradient search-^iethod x ^ x — j'S/H, 7 w 0.01. 

3) Ravine methodJla In this method the direction I = 
Xm — Xi of the most probable displacement of the system 
that moves in a long narrow ravine is predicted on the 
basis of a series {xi, X2, ■■■Xm} of previously performed 
gradient search steps. A number of trial moves in this 
direction x — + 0.5(a:i -t- X2) + IS completes the method. 
The maximal size of these moves max{|(5|} was adjusted 
in such a way as to achieve the acceptance probability of 
^0.1. We found that the throughoutput of the combined 
method was maximal when moves of above - mentioned 
types where performed in the ratio of Ni : N2 : N3 = 
100 : 10 : 2. 

When studying the thermodynamic proaerties of the 
system, we used the Metropolis algorithm^ For the sys- 
tem of small number of particles {N < 40) and at suf- 
ficiently low temperatures {T < 0.02), when shells are 
well defined, we found it very efficient to perform col- 
lective M|G, moves of different shells. With the multigrid 
approaches for shell s of Ns particles one can apply the 
following types of collective moves with wave vector ks : 

1) angular perturbations such that angles tpi^ of parti- 
cles belonging to shell s vary as 

ipz, Vis + £,S^{k; s) cos {2Trkis/Ns) (2) 

with ^ being the random variable uniform in [—1, 1). 

2) radial perturbations of shell s: 

Tis ^is + £,Sr{k] s) COS {2TTkis/Ns) (3) 



It is obvious that the case ks = corresponds to the 
symmetrical " breathing" of shell s when the radial per- 
turbation 2) is concidered and to the rotation of the shell 
as a whole when one performs the global move of the first 
type. We found that such rotation of different shells was 
very important when studying the phenomena of inter- 
shell orientational disordering which take place in meso- 
scopic systems at low temperatures (the temperatures 
of orientational "melting" can be many orders less than 
temperatures of full disordering and shells destruction). 
The use of the collective moves described above enabled 
us to increase an efficiency of calculations, that is in- 
versely proportional to the time needed to obtain a ther- 
modynamic average with a given accuracy, about a five 
times. 

Parameters and Sr (in Equations (H),(|^)), as well as 
the maximum size of the Metropolis trial move, were ad- 
justed in such a way as to hold the acceptance probability 
Pa « 0.4. To do this we averaged the acceptance prob- 
ability (p) of trial moves of different types {{pr) , {Pip)) 
in ~ 50 MC steps per particle followed by the scaling 
5^{k] s) 5^{k\ s)(l -Pa)l (p^) (and analogous for Sr). 

To study dynamical characteristics of the system we 
applied the molecular dynamic simulation (both isoki- 
netical and microcanonical methods where used). The 
main portion of the results presented below were obtained 
with the help of the fourth - order Runge - Kutta scheme. 
Equations of motion were integrat ed in ^ lO"' MD steps 
of size T < O.OItq, where tq = -y/ a/m* determines the 
time scale in the system. 



III. GROUND STATE CONFIGURATIONS 

studies of Coulomb00 and logarithmicii@0 clusters 
have shown that it is suitable to classify these finite sys- 
tems in accordance with their shell structure. An analy- 
sis of shell structures for different number of particles N 
enables one to consider the system as belonging to some 
period of a Mendeleev - type table. This table can be 
viewed as a classical equivalent to the well - known Peri- 
odic Table of elements. By studying regularities in this 
table the Coulomb clusters with a constant density of 
particles in each shell have been explored. This systems 
have been called "magic" Coulomb clusters. As exam- 
ples of such systems clusters of 18 (1, 6, 11), 30 (5, 10, 15) 
and of 57 (1,6,11,17,22) particles have been proposed. 
The existence of magic clusters particularly manifests it- 
self as cusps at a plot of excess energy €{N) = E{N)/N 
as a function of the number of particles. It have also 
been shown that these systems are more stable against 
intershell rotation. 

Interaction between particles in a dipole cluster is more 
short - ranged with respect to that in Coulomb or log- 
arithmic systems. Therefore it is possible to trace the 
role of the range of interparticle potential in the struc- 



3 



ture of 2D clusters, in the character of the formation of 
the Periodic Table and the type of shells (see also Ref. [ 
[Tot). ^ difference between structures of dipole (D) and 
Coulomb clusters can be seen even at = 10 (the sys- 
tem Dio has the structure 1510(3, 7) as distinct from (2, 8) 
for Coulomb system. o) The number of differences grows 
rapidly with increasing the number of particles. To ex- 
plore tendencies in the process of cluster shells develop- 
ing, let us consider the shell configurations of 2D dipole 
clusters that are presented in Table I. The table shows 
that the basis for most configurations is provided by dif- 
ferent parts of 2D hexagonal lattice. It should be noted 
that this observation is not specific to clusters with short 
- ranged interaction only. As was argued by Koulakov 
and Shklovskii [ |lo|, only a narrow ring adjacent to the 
perimeter of sufficiently large Coulomb clusters is concen- 
tric to it, the rest of the cluster is filled with an almost 
perfect crystal. 

When describing and analyzing the properties of such 
configurations we found it suitable to introduce into con- 
sideration the "crystal shells" Crj, that are concentric 
groups of nodes of ideal 2D crystal with c nodes placed 
in the center of these groups (see Fig. 1). Obviously, in 
view of the axial symmetry of the confinement poten- 
tial, we can concentrate on a finite number of the most 
symmetrical crystal shells. 

By the number of particles in the center of the sys- 
tem, the crystal shells can be divided into the following 
groups: Cri, Cr2, Cra, Cr4. Fig. 1 explains our defini- 
tion. The number of particles Ng that belong to crystal 
shell s (crystal row) of type Cvc is Ng = c + 6{s — 1) . 
With the help of the crystal shell concept an analysis of 
the results presented in Table I shows that a formation 
of a "Periodic Table of elements" is done in accordance 
with the following tendencies: 

1. The maximum number of crystal shells is filled. 

2. The number of particles on the last two shells tends 
to be equal. 

To illustrate these tendencies we underline the shells 
which are filled up (see Fig. 2 and Table I where types of 
basic crystal shells are also shown). 

An addition of particles to the cluster leads to the com- 
pleting of crystal shells of some type, followed by the 
structure rearrangement (see Fig. 2(c)-2(d)) after which a 
crystal shells of different type begins to fill. From Table I 
one can see that in some cases it is more advantageous to 
departure from the order Cri —f Cr2 Cr^ —^ Cr^ in 
which different types of crystal shells appear. Examples 
of such deviations (see D^q — Dgi and Dji ~ Dj^) show 
that it can be profitable to reduce the number of particles 
on first shells (in the center of the system) and equalize 
them on the last ones. This observation underlines the 
importance of the requirement that the last two shells of 
the cluster have equal number of particles. 

It is to be noted that for some systems the choice of the 
basis crystal group is ambiguous. The most spectacular 



examples of this ambiguity are clusters that we assign to 
Cr2 crystal group (such clusters are D^g, D^g, Dqi, Dq2, 
see Table I and Fig. 3). Of course, one can consider these 
systems as having a number of partially completed crys- 
tal shells of Cr4 type (the "body" of a cluster) while the 
other particles locate on the surface of this " body" . This 
approach is illustrated in Fig. 3 where two possible types 
of breakup of cluster into shells is shown. In choos- 
ing between these possibilities we have being guided by 
the requirement the outer shell of the cluster to be well 
- defined and closed. For the system with a more short - 
ranged interaction potential (e.g. with the exponential or 
screened Coulomb interaction potential) it can be more 
natural to consider the evolution of the cluster with N as 
consisting of switching between three types of the most 
symmetrical basic crystal shells, namely Cri, Cr^, Cr^&B 
Seemingly, the question of what approach is more ade- 
quate can be answered on the data of the thermodynamic 
and dynamic analysis. 

Results presented in Table I and above considera- 
tion enable us to suppose that for 2D dipole clusters 
a "period" of the Periodic Table consists of clusters 
with partially completed crystal shells of some symme- 
try (of some basic group). A role of "magic" clus- 
ters will be played by that with the maximal num- 
ber of completely filled crystal shells and, in the case 
of clusters of large number of particles {N > 40), 
with equal number of particles on the two last shells. 
An analysis of Table I suggests that dipole clusters 
Di2,Di4, Dig, Dse, D40, D^i^D^^, Dss, 1562, -Dgs... can be 
considered as "magic" ones. Plots of the first and the sec- 
ond "derivatives" of specific energy e{N) with respect to 
number of particles N (see Fig. 4) can serve as illustra- 
tions of this assumption. The particle numbers at which 
cusps take place correlate with the "magic" numbers. 

Somewhat ambiguous is the arrangements of particles 
into shells in the cluster 1)37 (see Fig. 5). We can see that 
one of the particles is between the second and the third 
shell to form an interstitial (it is analogous to the Frenkel 
defect in crystals). Such classification of this particle is 
based on the Voronoi analysis that shows that the cen- 
tral particle has six neighbors. The corresponding config- 
uration can be denoted by £'37(1, 6, 1, 13, 16) where the 
tilde stands for the interstitial. There exist also other op- 
portunities to describe the structure of this cluster, for 
example as 1)37(1,6, 12,2, 16). Our choice (see also Ta- 
ble I and Fig. 5(a)) is also proved itself in the results of 
the analysis of the lowest local minima configuration that 
has well - defined shell structure: d'"^^ {1, 7, 13, 16) This 
configuration is shown in Fig. 5(). It should be noted 
that the cluster involved has a rather complex picture of 
disordering phenomena with increasing temperature. In 
the study of its thermodynamic properties presented be- 
low we show that the system D^^ has properties that are 
characteristic of both small [N < 37) and large {N > 37) 
clusters. 

Another region of distinctions from the above tenden- 



4 



cies in particle distribution throughout shells is a group 
of clusters ZJyo — D75 that have a pentagonal shape. It is 
possible that the size AA'' of the region of non - hexag- 
onal arrangement of TV - particle cluster facets will be 
diminished with reducing the range of an interaction en- 
ergy. 

To obtain further information about the character of 
shells distribution in clusters of different number of par- 
ticles we applied the spectral analysis of the ground state 
configurations (the spectral analysis of Coulomb clusters 
was performed in the work [ ^). The main attention 
was given to the minimal nonzero eigenfrequencies ujmin 
for the normal modes. Fig. 6 shows oJmin as a function 
of the number of particles in the cluster. It is worthwhile 
to note that magic number clusters have maximal values 
of LOmin- The Correlation between the value of ojmm and 
the extent to which crystal shells are completed is also 
seen. 

One of the most interesting features of clusters is the 
strong dependence of their thermodynamic and dynamic 
properties on the number of particles N . An investi- 
gation of this dependence (see below) shows that dipole 
clusters can be divided into two groups: 1) "mesoscopic" 
("small") clusters, N < 37; 2) "macroscopic" ("large") 
clusters, N > 37. The intershell orientational disorder- 
ing phenomenon is distinctive of mesoscopic clusters and 
does not take place in dipole clusters of a large number of 
particles. But differences between large and small clus- 
ters can be already seen when considering ground state 
configurations (at T = 0). 

Before of all, increasing the number of particles leads 
to the changes in the structure of lowest excitations, i.e. 
in the picture of eigenvectors with smallest nonzero fre- 
quencies. Shown in Fig. 7 are some characteristic exam- 
ples of corresponding motions. Analogous to Coulomb 
clusters we found that small clusters with incommen- 
surate number of particles in neighbor shells have the 
smallest eigenfrequencies Wmm- The motion with mini- 
mal eigepfrequency LOmin is then corresponds to intershell 
rotationcil (see Fig. 7(a), (b)). When the number of parti- 
cles is increased one can see that the picture of a motion 
with the minimal eigenfrequency evolves in such a way as 
to fall into a number of individual regions of " rotations" 
to form "vortex/antivortex pairs" (Fig. 7(c), (d)). 

An interesting effect that accompany going from small 
to large clusters is the " facetting" . Shells of mesoscopic 
clusters have near - circular shape, while an addition of 
particles leads to the developing of well - defined facets. 
This tendency is clearly seen in Fig. 8 in which the de- 
pendence of the measure of facetting i/g on N for different 
cluster shells s is presented. If shell number s of the clus- 
ter consists of particles with numbers 12, --Iat}, then 
the parameter can be defined as 



max{|rjj} 

'is 

min{|rij} 



(4) 



of particles, the shape of the cluster shells approach that 
of the crystal shells defined above (appropriate values of 
the parameter Vs are also shown by thick solid lines) . In 
support of the fact that shells of small clusters have near 
- circular form. Fig. 8 shows that minimal differences be- 
tween values of parameter Vs for crystal and cluster shells 
take place when the most symmetrical crystal shell Cri 
is being filled (see Dn — Dig, D32 — D^q). Instead, this 
difference is maximal for clusters belonging to the other 
periods of " Mendeleev table" (i.e. ofCr2, Cra, Cr4 ba- 
sic groups). 



IV. PHASE TRANSITIONS 

Studying thermodynamic and dynamic properties of 
2D dipole clusters, we explored the temperature region 
T < 1 and calculated the following quantities: 

1) Radial distribution function g(r). 

2) Pair and radial deviations: 



i<J L 



1/2 



(5) 



1 ^ 



-,1/2 



(6) 



3) Radial diffusion constant Vr- This quantity changes 
drastically in the region of temperatures that is appro- 
priate to the radial disordering of the cluster and to the 
interchange of particles between different shells. For it 
is beyond reason to believe that a hydrodynamical ap- 
proach is valid for such small systems, we found it in- 
correct to define the diffusion constant via a Fourier - 
image of the velocity autocorrelation function. Instead, 
this quantity was estimated from equation 



(7) 



From Fig. 8 one can see that with increasing the number 



where (Ar^) = ^ Eili [k^ (01 " |r^(0)r is a radial 
mean square displacement of a particle in time t (in units 
of To, see above), c is some constant. 

4) Local minima distribution pieioc)- In order to es- 
timate this histogram and to determine the relative oc- 
cupancy of different minima on the potential energy sur- 
face, at each step of measurements several hundred steps 
of the gradient minimization (see above) was come out 
to determine the nearest local minimum. 

In simulations of macroscopic clusters we found it help- 
ful to calculate the following quantities: 

5) The number dn{Rc) of particles which are in a circle 
of radius i?c and have n nearest neighbors. The Voronoi 
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construction was applied to find nearest neighbors of each 
particle. 

6) Static structure factor S'(k) 

S(k) = ^ < PkP-k >, 

N 

Pk = Y^ exp (jkr,), f = -1 (8) 

i=l 

A usual way of studying the intershell orientational 
disordering in clusters is an analysis of relative angle de- 
viations of shells, in analogy with quantities (||) - (||). 
In this approach the temperature Ts-^S2 of orientational 
"melting" of shells si and S2 of a cluster is defined as 
that at which there is a sharp increase in the value of ap- 
propriate relative angle deviations. We follow a different 
method of finding the disordering temperature Ts-^S2 s-s 
that at which the mutual orientational ( order) parame- 
ter of shells si and S2 vanishes. We define this quantity 
as follows: for each shell s of particles we consider 
complex - valued quantity ij^s'- 

" ]^ 51 ^'^P (J^'*'^^) (9) 

^ is 

The sum in (||) is extended over all particles possessed 
by shell s. The mutual orientational (order) parameter 
is then defined as 

9s^s2 =< V'siV-:, > (10) 

It is obvious, that Qsis^ disappeares at the point of 
relative disordering (slipping) of shells si and S2 (in 
the vicinity of this the very point there will be also a 
sharp increase in relative angle deviations). The quan- 
tity gss —< IV'sP > will be a measure of the intrashell or- 
der. We note that quantities ips and gsiS2 a-re analogous 
to orientational parameter ipQ and correlation function 
ge{r) in infinite 2D system, where vanishing (when the 
translational order is absent) of the correlation function, 
56(^) ^ as r — !■ oo, indicates the relative orientational 
disordering of distant parts of a system. 

A. Mesoscopic clusters 

The distinctive property of mesoscopic clusters is the 
presencQ— p[L|k«o types of disordering effects in these 
systemsO'ESE^ an intershell orientational disordering 
(an orientational melting of shells si and §2 at temper- 
ature Tsjsj) and a radial disordering (a total melting at 
temperature Tf that is larger than any one of orienta- 
tional melting temperatures). An analysis of eigenfre- 
quencies and eigenvectors for the normal modes of small 
clusters shows that the motions with small lowest nonzero 
eigenfrequencies ujmin correspond to intershell rotation 
(see Fig. 7). Such clusters will have small temperatures 
TsiS2 of intershell disordering, at which shells start to 



rotate relative to each other losing their mutual orienta- 
tional order. Note that, in contrast to the case of large 
clusters, an intershell melting in small clusters takes place 
for all pairs of shells, i.e. there exist "melting" temper- 
atures T21, T^2i 743... . In clusters of large number of 
particles {N > 37) an orientational disordering of only 
the outer shell is possible. Aa. analogous observation was 
made for Coulomb clusters. EJ 

Shown in Fig. 9(b) are dependencies of the values of 
mutual orientational parameters ^21 and 332 vs. tem- 
perature for three - shell cluster i?24- Pair and radial 
deviations 5pair{T) and u^(T) (see (|1),(^)) are also plot- 
ted. This figure shows that for cluster D24 one can 
define two temperatures r32 = (5 ± 0.5) • 10^** and 
T21 — (3±0.05)-10~^ that are correspond to orientational 
melting of shells {3, 2} and {2, 1}. 

Fig. 9(b) (see also Fig. 10) shows that pair and radial 
deviations (|5|), (^ are also sensitive to intershell rota- 
tion. One can see the regions of sharp increases in the 
values of these quantities that coincide with the regions 
of vanishing of mutual orientational order. The sensi- 
tivity of radial and pair deviations to the orientational 
disordering is due to the " breathing" of the cluster shells 
on their rotation. This breathing can be clearly seen 
if onejirace the motion of a system along a "reaction 
path" ,c3 the most probable (of the weakest ascent in 2A^ 
- dimensional space) trajectory which is appropriate to 
the intershell rotation. 

Temperature T32 ~ 0.05 ± 0.003 of the orienta- 
tional "melting" of the third shell of four - shell cluster 
1535(1,6,12,16) is only slightly lower than the temper- 
ature Tf = 0.065 ± 0.005 of the total melting at which 
the radial order disappears (see Fig. 10). It is obvious 
that the temperature of an orientational intershell disor- 
dering will to a large extent depend on the distribution 
of particles throughout shells. Particularly, when the 
particle numbers in neighbor shells si and S2 are com- 
mensurate, temperature Ts^s2 of an orientational "melt- 
ing" should be maximal (irrespective of whether these 
shells are inside the cluster or a pair of external shells 
is considered). Calculations do show that an addition 
of as little as one particle to small cluster Dig(l,6,12) 
with completely filled crystal shells of Cri type leads 
to abrupt decreasing of temperature T32 of orientational 
"melting". Indeed, temperature T32 « 0.03 ± 0.002 of 
Dig cluster is very nearly equal to that of the total melt- 
ing Tf w 0.038 ± 0.003 (at which particles begin to in- 
terchange between shells). System -D2o(l, 7, 12) becomes 
orientationally disordered at much smaller temperatures 
r32 < 5 • 



B. Macroscopic clusters 

The most interesting when considering macroscopic 
clusters (at N > 37) is the question about the manner in 
which their melting temperatures Tf approach the tem- 
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peratureT™^ = kbT/{(fn^^'^) = 0.089 ± 0.002 of first or- 
der phase transition in 2D infinite dipole system |(bcre n 
stands for the density of particles in the system) UVl But 
this analysis is made difficult by an uncertainty in the 
concept of the "density" of a cluster. This uncertainty 
does not enable us to use the natural for infinite systems 
unit of length a — 1/ yfn. When converting from " clus- 
ter" to "infinite system" temperature, one can introduce 
a mean density of a system as n = N/-kR^ where R is the 
radius of the cluster. The uncertainty in R leads to the 
ambiguity of such definition of the cluster density. More- 
over, if we deal with a cluster of a large enough number 
of particles the density will considerably vary with the 
distance to the center of the system. We will discuss this 
question in some more detail below. 

To define an effective (mean) density of a cluster of 
not very large number of particles {N < 50), when the 
nonuniformity of the density pointed out above is not es- 
sential, we made use of calculations of the static structure 
factor (||). Let us consider the distance |<5k| to the maxi- 
mum of the structure factor, i.e. to the first Bragg peak 
for a "solid" cluster (when T < T/) or to the Lorentz 
peak in the case of a "liquid" cluster (when T > Tf). 
The characteristic distance between particles in a clus- 

1/2 

ter can be defined then as a = 27r/|(5k| — [2/(nV3)] 
In such a manner the corresponding temperature of 2D 
infinite system can be estimated as 

r^nf _ hT _ 1 \\Sk\Y . . 

£)2„3/2 ~ 0.8 [ 27r J ^ ^ 

It turned out that for systems D^q — D^2 the estimation 
T « r™/ works with a reasonable accuracy. 

One of the most representative quantities for 2D infi- 
nite dipole system is excess energy e that exhibits a jumg. 
of Ae w 0.04 at a temperature of phase transition T™'^.Q 
Our calculations show (see Fig. 11) that such a sharp in- 
crease in an excess energy does not take place at least for 
clusters with N < 50. 

Yet another quantity that is commonly used in ana- 
lyzing phase transitions in infinite systems is structure 
factor S'(k), Eq. (||). The lattice to liquid transition is 
identified through the vanishing of the first Bragg peak 
S'(k), k f« qi. We found that the magnitude of this 
peak ( the maximum of the structure factor in the region 
|k| > fcc ~ tt) is acutely sensitive to the disordering in 
clusters. The peak at small wave vectors |k| < tt that 
exists at any temperature due to finite size of the system 
involved is of no interest. Some characteristic examples 
of the behavior of the structure factor as a function of 
temperature are given in Fig. 12. 

An independent quantity a sharp change in which testi- 
fies about the order/disorder transition is the radial dif- 
fusion constant Vr, see equation (|^). Fig. 13 presents 
the results of calculations of this quantity, performed for 
clusters I?3g , DiQ and . Each data point in this figure 
was obtained by averaging over results of 5 independent 
experiments in which the evolution of the system was 



observed in a time that is typical for the diffusion of a 
particle for a distance of Ar ^ 1. The result of one of 
such experiments is shown in the insert of the figure. 

An analysis of temperature dependencies of above 
quantities (^ - (H) enables us to plot the total melt- 
ing temperature Tf{N) as a function of the number of 
particles in the cluster. This curve is shown in the in- 
sert of Fig. 14. Also shown are temperatures Tj''^ [N) 
of corresponding infinite system that are estimated with 
the help of Formula ([ll|). Fig. 14 enables us to point 
out some characteristic features of the N dependence of 
the disordering temperature, i) Melting of all clusters 
studied takes place at temperatures much more smaller 
than temperature T™-^ « 0.1 of the phase transition in 
an infinite system, ii) The most interesting peculiarity of 
function Tf{N) is that it is nonmonotonic. Most likely, 
this peculiarity is primarily connected with the fact that 
a cluster of N particles can be considered as a part of a 
strongly deformed crystal lattice. Hence one studies the 
melting of a part of imperfect crystal with "frozen- in" 
interstitials and dislocations (see Fig. 2(a)), the melting 
temperature being a function of the number of such de- 
fects (of the value of an initial deformation of a cluster 
as a part of crystal). Obviously, the number of defects in 
the ground state configuration will be minimal in the case 
of " magic" cluster (see above) with a maximal number of 
filled crystal shells. An analysis of the data presented in 
Table I and in Fig. 14 shows that (at least for N < 52) 
the "magic" clusters do have anomalously high melting 
temperatures. 

Of course, the melting temperature is not a strictly de- 
fined quantity for finite systems. It makes sense to con- 
sider an interval AT{N) of temperatures in which the 
disordering takes place. For clusters of large number of 
particles that repel each other and are in a confinement 
potential there is another, more important, reason that 
does not enable us to introduce a concrete melting tem- 
perature of a phase transition (of a total disordering). 
The point is that large clusters (at N > 50) are irregular 
in both the density of particles and the density of defects. 
These densities can vary appreciably with the distance r 
to the center of the systemE3 to lead to a strong r depen- 
dence of the melting temperature. 

Let us consider this dependencies in some more detail. 
Results of numerical simulations show that large clus- 
ters can be divided into the following regions £3 (a) the 
"core" of a cluster, a region adjacent to the center of 
the island and filled with a number of completed crys- 
tal shells; (b) the "layer of defects", the ring in which 
dislocations and disclinations are situated and which is 
"concentric" to the surface of the cluster; (c) the re- 
gion between the surface of the island and the layer of 
defects. One can assume that temperatures of disorder- 
ing in the regions considered above will differ greatly. 
Our calculations approve this assumption. As an exam- 
ple, we calculated the number of disclinations (i„(r) and 
the radial diffusion constant Vrir) as functions of the 
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distance r to the center of the cluster Dso- It turned 
out that temperature Tj^^ {f < ^c) of the disordering 
of a "core", the part of the cluster belonging to the 
circle of the radius Vc ~ 2.5, was independent of r: 
Tp-^{r < 2.5) « 0.07±0.01. This temperature is slightly 
below T^"^ . One may suppose that as N increases, the 
scaled temperature of core disordering Tp-^{N;r < Vc) 
approaches the temperature of phase transition in 2D in- 



N ■ 



oo. 



finite dipole system: Tp {N; r < r^) ^ T^''^ , 

The melting in other two regions of the cluster ap- 
peared to take place at much less temperature: T™^ (r > 
Tc) = 0.01 ± 0.005. Obviously, the temperature of dis- 
ordering in regions (b) and (c) will strongly depend on 
the number of disclinations and dislocations in the layer 
of defects (in the region (b)) as well as on the strength 
of initial deformations in the region (c) . Thus we argue 
that: 

i) defined as a point of drastic changes in pair and ra- 
dial deviations, diffusion constant and of disappearance 
of the first Bragg peak, neither "melting" temperature 
Tf (N) nor scaled to 2D infinite system melting tempera- 
ture Tj''^ {N) tends to first - order phase transition tem- 
perature in infinite 2D system as the number of particles 
in the cluster increases. 

ii) The function Tf^{N) (as well as Tf{N)) is not mono- 
tonic. The magic number clusters, having minimal num- 
ber of defects, have maximal melting temperatures. 

iii) One can introduce the temperature Tj^^{N;r < 

Tc) T^"-^, N ^ oo which is the (scaled to 2D in- 
finite system) melting temperature of a core, the free of 
defects central region of a cluster. 
These qualitative predictions are presented in Fig. 14. 



C. Cluster D37 

As have been pointed out above, system D37 can be 
considered as an intermediate one between mesoscopic 
and macroscopic clusters. An analysis of the global min- 
imum configuration (see Fig. 5) and of the structure of 
normal motion with the minimal eigenfrequency ujmin 
leads to the seemingly unambiguous conclusion about the 
absence of the intershell orientational disordering. How- 
ever, our calculation shows that the melting in this sys- 
tem is not one - step process; instead, the orientational 
melting takes place at the temperature T = T* k, 0.01 
that is approximately three times lower than that of total 
melting Tf m 0.03 ± 0.005. The orientational disordering 
of all pairs of shells takes place at temperature T* . 

The reason of this "anomalous" behavior of cluster 
1^37 can be elucidated by considering the structure of 
the lowest local minima configurations (i.e. the config- 
urations with the lowest energy). We have noted above 
(see Fig. 5(c)) that global and local minima configura- 
tions differ in a symmetry^ namely, any breaking of the 
ground state configuration down into shells leads to the 



necessity of viewing one or two particles (see Chapter III) 
as interstitials interposed between shells of the cluster. 
Only two first shells are well defined: 1)37(1,6,...). On 
the contrary, the configuration of the first excited state 
(of the lowest local minima with energy e|^^ ~ 8.3325) 
is rather symmetrical and has the clear shell structure: 
Dgy-* (1, 7, 13, 16). As shells of the cluster in this local 
minimum are not facetted, it is not surprising that the 
motion (in this local minimum) with the lowest eigenfre- 
quency i^£]„ ~ 0.4 corresponds to the intershell rotation. 
This is illustrated in Fig. 5(c). 

Plotted in Fig. 15 is the probability of the central par- 
ticle having six neighbors. This quantity can be defined 
as (i6(0.5), see above. Fig. 15 shows that at tempera- 
tures T >T* = 0.01 ± 0.002 there is a finite probability 
to find a system in the vicinity of the local minimum 
£'37-' (1, 7, ...). Yet another illustration of this rearrange- 
ment is the change in histogram p{eioc) of the local min- 
ima distribution shown in the insert of Fig. 15. One can 
see that at T = T* a pike at energy e|^' appears. Fur- 
ther increasing the temperature leads to the occupancy 
of others local minima and, at T > 0.3, the radial dis- 
ordering takes place, particles interchange between shells 
and the radial distribution function is washed out (see 
Fig. 15,16). So, at r > T*, a part of time the system 
lives in the neighborhood of the " symmetrical" local min- 
imum that is unstable against intershell rotation at such 
high temperatures. This intershell disordering also mani- 
fests itself as a sharp increase in relative angular intershell 
deviations at this the very temperature T = 0. 01 . ■ 

From results of studies of atomic clusterai^liil it is 
known that clusters of some specific number of particles 
(so called "magic" LJ clusters such as LJ13, iJig, see 
work [ ^] and references therein) have distinct regions of 
coexistence of " liquidlike" and " solidlike" forms. In a cer- 
tain temperature interval there is nonzero probability to 
find the cluster either in a " liquidlike" or in a " solidlike" 
form. It was shown that this unusual property owes its 
origin to features of the multidimensional potential en- 
ergy surface and stems from an existence of regions (well 
- separated by high barriers) in phase space with distinct 
physical properties and appreciable residence time for the 
system in each region. 

A picture of disordering in cluster D^^ with increasing 
temperature closely resembles the picture of structural 
transitions in magic number atomic clusters. Of course, 
at sufficiently high temperatures, at T > T* = 0.01, the 
system can be found in one of two distinct regions in 
phase space with a barrier between them. This equilib- 
rium can be described as: 



K 



i?37(l,6,l,13,16)^i?37(l,7,13,16). 



K 



[i^37(l,7,13,16)] 



^-AF{T)/T 



[i?37(M'l'13,16)] 
where AF{T) is free energy difference of "symmetrical" 
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and "fully ordered" states of the cluster. At T < 0.01 
constant of equilibrium K ^ \ — ^5(0. 5) is equal to zero 
and increases markedly as the temperature is increased 
(see Fig. 15). 

The analysis above shows that, with the cluster D^-j 
being radial ordered, the region 0.01 < T < 0.03 can not 
be considered as the region of coexistence of its "solid- 
like" and "liquidlike" forms. Instead, this temperature 
interval can be though of as the region of coexistence of 
"solidlike" and "orientationally disordered" forms. It is 
the presence of such region that enables us to consider 
the system D37 as an intermediate one that has features 
of both mesoscopic and macroscopic clusters. 

To conclude of this paragraph we note that the number 
of 1)37 - like " anomalous" systems can be varied by tun- 
ing the characteristic range of an interaction potential. 
Increasing the temperature in such systems can lead to 
structure rearrangements (with changes in symmetry), 
changes in the shell distribution or in types of defects. 
Consequently, above consideration can be helpful in the 
study of thermodynamic properties of the clusters, pa- 
rameters of interparticle potential of which can her-war- 
ied in wide limits, e.g. of dusty plasma clusters, cl-Ea of 
the system of electeens in a semiconductor dot with near 
metal electrodes,Q'liy etc. 



The system of large number of particles is nonuniform: 
both the characteristic interparticle distance and the den- 
sity of defects are functions of a distance r from the cen- 
ter of the system. As a consequence, the temperature 
at which a sharp increase in radial and pair deviations, 
vanishing of the first Bragg peak take place and the dif- 
fusion of particles appear does not approach temperature 
of first - order phase transition in 2D infinite dipole 
system as the number of particles is increased. Scaled 
to infinite system, "melting" temperature Tj^^ {N) is a 
nonmonotonous function of the number of particles and 
is maximal for magic number clusters. By taking into ac- 
count the r dependence of the interparticle distance, the 
melting of the "core" of the cluster (i.e. its free of defects 
central region r < rdN)) can be mapped to the melting 
of 2D infinite system with appropriate scaled tempera- 
ture T'^'fiN; r < Tc) -> T^'^^ as N increases. 
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V. CONCLUSION 



We have presented the results of a numerical simula- 
tion of finite 2D dipole system in the parabolic confine- 
ment. Ground state configurations and the spectrum of 
normal modes of clusters of iV < 80 particles have been 
found. An addition of particles to the system leads to a 
gradual filling of crystal shells of different symmetry. The 
clusters with minimal energies and maximal lowest eigen- 
frcquencies (" magic" clusters) are that with the maximal 
number of completely filled crystal shells and (in the case 
of large clusters, N > 40) with equal numbers of parti- 
cles on the last two shells. Such clusters have also well 
developed facets. 

The character of disordering with increasing the tem- 
perature is markedly different for mesoscopic {N < 37) 
and macroscopic {N > 37) clusters. Particularly, meso- 
scopic clusters are characterized by an existence of two 
types of disordering effects: orientational and full (ra- 
dial) . Depending on the degree to which crystal shells are 
completed the temperatures of orientational " melting" of 
different pairs of shells can differ greatly from each other. 
Orientational melting in large systems is absent. 

An analysis of the local minima distribution of the sys- 
tem Z?37 as a function of temperature have shown that 
this system has features of both mesoscopic and macro- 
scopic clusters. Namely, there is a region of temperatures 
in which cluster is in a dynamical equilibrium between 
"solidlike" and radial ordered but orientationally disor- 
dered forms. 
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Table I. 

Ground state configurations of 2D dipole clusters Djv in a harmonic confining. 
{Ni,N2, ...}, types of crystal shells (see in Text) and excess energy e = E/N. 



Shown are shell configurations 



N 


{Ni,N2,...} 




e 


N 


{Ni,N2,...} 




e 


1 


1 







41 


3,9,14,15 


Crs 


8.93397 


2 


2 




0.64660 


42 


3,9,14,16 


Crs 


9.08148 


3 


3 




1.01394 


43 


3,9,15,16 


Crs 


9.22705 


4 


4 




1.38021 


44 


3,9,15,17 


Crs 


9.37309 


5 


5 




1.75713 


45 


4,10,15,16 


Cn 


9.51691 


6 


1,5 




2.04829 


46 


4,10,15,17 


Cn 


9.66006 


7 


M 


Cri 


2.32591 


47 


4,10,16,17 


Cri 


9.80303 


8 


1,7 


Cri 


2.63542 


48 


4,10,16,18 


Cn 


9.94353 


9 


2,7 


Cr2 


2.92373 


49 


1,5,11,16,16 


Cri 


10.08318 


10 


3,7 


Crs 


3.19012 


50 


1,6,11,16,16 


Cn 


10.21803 


11 


3,8 


Cr:i 


3.41972 


51 


1,6,12,16,16 


Cri 


10.35593 


12 


3,9 


Cra 


3.66665 


52 


1,6,11,17,17 


Cn 


10.49035 


13 


4,9 


Cr-4 


3.89493 


53 


1,6,12,17,17 


Cn 


10.62096 


14 


4,10 


Cr-4 


4.13543 


54 


1,6,12,17,18 


Cn 


10.75525 


15 


5,10 




4.35999 


55 


1,6,12,18,18 


Cri 


10.88617 


16 


1,5,10 




4.56558 


56 


1,6,12,18,19 


Cn 


11.02094 


17 


1,6,10 


Cri 


4.77272 


57 


1,6,12,18,20 


Cn 


11.15521 


18 


1^,11 


Cri 


4.97257 


58 


1,6,12,18,21 


Cn 


11.28939 


19 


1,6,12 


Cri 


5.18009 


59 


2,8,13,18,18 


Cr2 


11.41878 


20 


1,7,12 


Cn 


5.38833 


60 


3,8,13,18,18 


Crs 


11.54741 


21 


2,7,12 


Cr-2 


5.59048 


61 


2,8,14,18,19 


Cr-2 


11.67875 


22 


2^,12 


Cr2 


5.77969 


62 


2,8,14,19,19 


Cr-2 


11.80329 


23 


3,8,12 


Crg 


5.96866 


63 


3,8,14,19,19 


Crs 


11.92863 


24 


3,8,13 


Crs 


6.14713 


64 


3,9,14,19,19 


Crs 


12.05160 


25 


3,9,13 


Crs 


6.32561 


65 


3,9,15,19,19 


Crs 


12.17598 


26 


4,9,13 


Cr4 


6.50834 


66 


3,9,14,20,20 


Cr-s 


12.30108 


27 


4,9,14 


Cn 


6.68410 


67 


3,9,15,20,20 


Crs 


12.42251 


28 


4,10,14 


Cr4 


6.85654 


68 


3,9,15,20,21 


Crs 


12.54674 
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Fig. 1 

Different parts of 2D hexagonal lattice produce the basis 
for ground state configurations of the majority of dipole 
clusters (see Table I). Shown are examples of such basic 
groups of crystal shells. Four most symmetrical groups 
of three shells each is presented. 

Fig. 2 

Ground state configurations of dipole clusters: a) -D44; 
b) Z?47; c) Z?52; d) Z^sg. Squares denotes particles with 5 
nearest neighbors. 

Fig. 3 

Two possible breakups of cluster 0^2 into shells: a) in 
such a-syay as to maximize the number of completed Cr^ 
shells.E3 b) so that the outer shell be well - defined and 
closed. 

Fig. 4 

The second and the first (in the insert of the Figure) 
"derivatives" of specific energy e{N) with respect to N. 

Fig. 5 

Cluster I?37. a) The ground state configuration 
1)37(1,6, 1, 13, 16). b) the eigenvector of the cluster mo- 
tion in the global minimum with the minimal nonzero 
eigenfrequency Wmm ~ 0.58 c) the picture of the motion 
in the lowest local minimum 13^7^(1,7,13,16) with a min- 
imal eigenfrequency « 0.4. 
Five - coordinated particles are marked by squares. 

Fig. 6 

Minimal nonzero eigenfrequencies for the normal 

modes of - particle dipole clusters. 

Fig. 7 

Characteristic examples of normal motions with eigen- 
frequencies uj,nin (see Fig. 6). a) Diq; b) D35; c) D55; d) 

^56- 

Fig. 8 

The measure of faceting i's(iV) is shown for different 
shells s. Appropriate values for crystal shells of differ- 
ent types are plotted by thick solid lines. 

Fig. 9 

Three - shell cluster £)24(3, 8, 13). The results of calcu- 
lations of the mutual otientational parameter ( p^ and 
of deviations (||),(^) as functions of temperature T are 
presented. 

The data are connected to guide the eyes. If not present, 



error bars are smaller than the size of the data point. 
Fig. 10 

Four - shell cluster I?35(l, 6, 12, 16). Mutual orientational 
parameter _Re|(732| and relative deviations 6pair(T) vs. 
temperature T. 

Fig. 11 

Temperature dependence of excess energy e(T) for a num- 
ber of dipole clusters Dn. In the insert the excess energy 
of D37 cluster in the temperature region T E [0.02,0.02] 
is shown. 

Fig. 12 

The magnitude of the first Bragg peak as a function of T 
for clusters £'40 and -D45. Shown also is 2D topography 
of the structure factor S'(k), S'(k) < 0.5 at different 
temperatures. "Solidlike" to " liquidlike" transition is ac- 
complished by the washing out of S(k) in regions k « qi 
of reciprocal lattice vectors. 

Fig. 13 

The radial diffusion constant (0) vs. temperature T. The 
insert shows a typical diffusive motion of particles when 
the cluster is in disordered state (T > Tf). Results of 
the linear fit (see Equation ^) are also given. 

Fig. 14 

The total melting temperatures Tf{N) and scaled from 
" cluster" to " infinite system" temperatures Tj^-^ {N) ( |Tl| ) 
as functions of number of particles N. Points correspond 
to the results of present simulations. The dotted line 
depicts a possible behavior of {N) a,t N > 50. Tem- 
perature of first - order phase transition in an infinite 2D 
dipole system T^""^ is shown with the help of a dashed 
line. Temperature Tj''^[N]r < Tc) of the disordering of 
a core is assumed to tend to T™-^ as N increases. 

Fig. 15 

The probability that the central particle of cluster D37 
has six neighbors as a function of temperature T . In the 
insert local minima distributions for different tempera- 
tures are shown. One can see that at T > 0.01 a group 
of local minima with energies near e^^-* , the energy of the 
lowest local minimum, are being occuped. 

Fig. 16 

Radial distribution function g{r) of cluster D37. 
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